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Background: The rate App^ characterizes the formation of ppp, molecules in collisions of muonic pfj, atoms with 
hydrogen. In measurements of the basic weak muon capture reaction on the proton to determine the pseudoscalar 
coupling gp, capture occurs from both atomic and molecular states. Thus knowledge of App^ is required for a 
correct interpretation of these experiments. 

Purpose: Recently the MuCap experiment has measured the capture rate As from the singlet pp atom, em¬ 
ploying a low density active target to suppress ppp formation (PRL 110, 12504 (2013)). Nevertheless, given the 
unprecedented precision of this experiment, the existing experimental knowledge in App^ had to be improved. 

Method: The MuCap experiment derived the weak capture rate from the muon disappearance rate in ultra-pure 
hydrogen. By doping the hydrogen with 20 ppm of argon, a competing process to ppp formation was introduced, 
which allowed the extraction of Appp from the observed time distribution of decay electrons. 

Results: The ppp formation rate was measured as Appp = (2.01 ± O.OOstat ± 0.03ays) x 10® s“'^. This result 
updates the App^ value used in the above mentioned MuCap publication. 

Conclusions: The 2.5 x higher precision compared to earlier experiments and the fact that the measurement 
was performed at nearly identical conditions to the main data taking, reduces the uncertainty induced by Appp 
to a minor contribution to the overall uncertainty of A 3 and gp, as determined in MuCap. Our final value for 
Appp shifts Ag and gp by less than one tenth of their respective uncertainties compared to our results published 
earlier. 

PACS numbers: 23.40.-s, 13.60.-r, 14.20.Dh, 24.80.+y, 29.40.Gx, 36.10.-k 


I. INTRODUCTION 


Nuclear muon capture on the proton, 

p~ + p^n + Uf_, , ( 1 ) 

is a basic charged-current weak reaction Several 

experiments have measured the rate of ordinary muon 
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capture (Eq. 0 ) or the rarer process of radiative muon 
capture, p+p —> n+i/+j, in order to determine the weak 
pseudoscalar coupling of the proton, gp, which can be ex¬ 
tracted most straightforwardly from muon capture on the 
nucleon. A precision determination of gp has been a long¬ 
standing experimental challenge Elia due to the small 
rate of capture on the proton and complications arising 
from the formation of muonic molecules. The most recent 
MuCap result, gp = 8.06 ±0.55 [ 1 ], achieved an unprece¬ 
dented precision of 7 %, thereby providing a sensitive test 
of QCD symmetries and confirming a fundamental pre¬ 
diction of chiral perturbation theory, gp = 8.26±0.23 [3- 

13 . 

Experimentally, process ([T]) is observed after low- 
energy muons are stopped in hydrogen, where they form 
pp atoms and ppp molecules. The overlap in the wave- 
functions of the proton and the bound muon leads to 
small but observable capture rates at the 10“^ level rel- 
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ative to muon decay, fi~ —> e~VeV^, which is the domi¬ 
nant mode of muon disappearance in that environment. 
The nuclear capture rates depend on the spin composi¬ 
tions of the muonic atoms and molecules (a direct con¬ 
sequence of the V—A structure of the electroweak inter¬ 
action), and thus the rates vary significantly among the 
different muonic states. The calculated rates for the two 
hyperfine states of the p/r atom possessing spin F=0,1 
are Ag = 712.7 s“^ and Arp = 12.0 s“^, respectively 
(c.f. [5, updated in m)- The formation of ppfi molecules 
further complicates the situation, as the calculated cap¬ 
ture rates for the ortho and para states, Aq,^= 542.4 s“^ 
and A_pm= 213.9 s“^, differ too from the atomic rates 
(Eq. Q). Correct interpretation of the observed muon 
disappearance rate thus relies on a thorough understand¬ 
ing of the “muon chemistry” reactions governing the time 
evolution of the pp and ppp states. This interrelationship 
between muon capture and muon chemistry in hydrogen 
has been the primary source of ambiguity in the 50-year 
history of experiments in the field. Historically, inter¬ 
est in muon atomic and molecular reactions arose due 
to their above mentioned relevance for the determination 
of nuclear muon capture rates in hydrogen isotopes [5] 
and their importance in muon-catalyzed fusion [5] , where 
Xpp^ was calculated within a systematic program to solve 
the Coulomb three-body problem [9]. 

The MuCap experiment employed a novel technique 
involving the use of low-pressure hydrogen gas to sup¬ 
press molecular formation. Nevertheless, it was still nec¬ 
essary to apply corrections which were based on measure¬ 
ments of the molecular formation rates that determine 
the ppp ortho and para molecule populations. In the 
initial MuCap physics result m, we conservatively es¬ 
timated that the uncertainty in the molecular formation 
rate contributed a systematic uncertainty of 4.3 s“^ 
to our determination of Ag, the muon capture rate in the 
pp hyperfine singlet state. During the later high statistics 
data taking for MuCap, we performed a dedicated mea¬ 
surement of Appp in order to improve the precision on 
this parameter and render its contribution to the uncer¬ 
tainty on Ag nearly negligible. The final MuCap result, 
Ag = (714.9 ± 5.4stat ± S.lsyst) gj, possessed greatly 
improved statistical and systematic uncertainties. A pre¬ 
liminary value for Xpp^ obtained from our measurement 
was an important ingredient to this result. In this paper 
we document the Xpp^ experiment and present its final 
results. 


The contents of this article are as follows: In Section HIl 
we introduce muon-induced processes in hydrogen and 
their impact on muon-capture measurements. In Sec- 
tion |III| we describe the MuCap experiment and our tech¬ 
nique for measuring App^; the corresponding data anal¬ 


ysis and result for App^ are described in Section IV In 
Section lv| we use our new result for App^ to update pre- 
vious AluCap measurements. A concluding summary is 
given in Section |vg 


II. MUON CAPTURE AND MUON 
CHEMISTRY 


A. Muon Reactions in Hydrogen 


Muons stopped in hydrogen can form a variety of 
atomic and molecular states which are subject to dif¬ 
ferent physical processes and whose populations are gov¬ 
erned by the rates shown in Fig. [2 Table 1^ lists all of the 
rates used in this paper and their values. Several of the 
atomic processes proceed via binary collisions of muonic 
atoms with other target molecules. It is conventional to 
normalize those density-dependent rates to the values ob¬ 
served at LH 2 density, (f)Q = 4.25 x 10^^ atoms/cm^, and 
express all target densities (f) relative to ^q. 



FIG. 1. (Color online) Kinetics of negative muons in hydrogen 
containing a single chemical impurity Z. The large circles rep¬ 
resent the different muonic states that can form. The black 
arrows denote transitions between muonic states, while the 
colored arrows with a small circle at the beginning of the line 
indicate muon disappearance due to the weak-interaction pro¬ 
cesses of decay or nuclear capture. The arrow labeled Apf is 
dashed to indicate that its rate is about 240 times smaller 
than Aof. 


Negatively charged, low-energy muons entering hydro¬ 
gen are slowed down and undergo atomic capture, form¬ 
ing highly excited pp atoms. After an atomic cascade 
to the ground state, the two hyperfine states of the pp 
atom, singlet (F=0) and triplet (F=l), are populated ac¬ 
cording to their statistical weights \ and |, respectively. 
These complex initial stages happen on a timescale of 
nanoseconds at target densities exceeding </> > 0.01, as in 
our case. Charge-exchange collisions EBUa convert the 
higher-lying triplet state to the lower-lying singlet state 
at a rate calculated to be « </> x 2 • 10^° s“^ [^. Thus 
after less than 100 ns the triplet state is effectively de¬ 
populated and the main features of the kinetics can be 
described by the scheme depicted in Fig. This condi¬ 
tion is true for the present analysis and in Refs. min]- 
For our purposes, muon kinetics in pure hydrogen ef¬ 
fectively starts with the pp atom in its hyperfine sin- 
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TABLE I. Summary of relevant rates for processes involving 
muons in hydrogen. 


Muon process 

Symbol 

Value [s ®] 

Reference 

Weak-interaction rates 




Muon decay 

A+ 

455170.05 ±0.46 

[nHS] 

pp singlet capture Ag 

714.9 ±7.4 

m 


pp triplet capture A^ 

12.0 ±0.1 

[D 


ppp ortho capture Aqm 

544.0 ± 11.3 

Eq. 


ppp para capture Api^j 

214.6 ±4.2 

Eq. 


N capture 

An 

6.93 ±0.08 X 10® 

m 


0 capture 

-^0 

10.26 ±0.06 X 10® 

m 


Ar capture 

AAr 

141 ± 11 X 10® 

m 




130.2 ±3.2 X 10® 

this work 

Atomic and molecular rates 




o-p transition 

"^00 

6.6 ±3.4 X 10® 

m 


ppp formatior|‘‘ App^ 

2.3 ±0.2 X 10® 

m 




2.01 ± 0.07 X 10® 

this work 

Transfer to h ® 

ApN 

0.34 ±0.07 X 10®® 

m 


Transfer to C ^ 

Apo 

0.85 ±0.02 X 10®® 

m 


Transfer to Ai 

d ApAr 

1.63 ±0.09 X 10®® 

m 



1.94 ±0.11 X 10®® 

this work 

Ar Huff factoi 

® h 

0.985 ±0.003 [IBJinj 


^ normalized to LH 2 density (f>Q 
dimensionless quantity 


glet state. In subsequent collisions of the pfi atom with 
hydrogen molecules, two types of ppp molecules can be 
formed which differ in their angular momentum L and to¬ 
tal spin I. Due to the Fermi statistics of the two-proton 
system, the ortho state {ppp)^^ has L=l, 1=1, while the 
para state {ppp)^^ has L=0, 1=0. According to theory, 
ppp, formation proceeds to the ortho state predominantly 
at the normalized rate = 1.8 x 10® s“^, while the para 
formation rate Apf=7.5 x 10® s“® is much smaller [5]. The 
total normalized molecular formation rate is the sum of 
these two rates, 

^ppfl '^of ^pf ■ (^) 

Molecular formation scales with the target density </>, 
so experiments observe the effective molecular formation 
rate 

\pp, = (l>\pp ■ (3) 

The transition from the ppp ortho state to the lower 
para state at rate A^p involves a proton spin flip and 
is only allowed due to relativistic effects in the molecu¬ 
lar wave function. The (pppom)^ is positively charged 
and quickly forms various molecular complexes in colli¬ 
sions with H 2 molecules. The ortho-para transition pro¬ 
ceeds at the calculated rate A^p = (7.1 ± 1.2) x 10® s“^ 
via the emission of an electron from these clusters 
Two previous experiments measured Agp and obtained 
the inconsistent results (4.1 ± 1.4) x 10® s“® [53] and 
(ll.l±1.9)xl0®s-® [ig. Review [2] therefore inflated the 
uncertainties and quoted an average experimental value 
of Aop = (6.6 ± 3.4) X 10® s“®, which we use in this work. 


As mentioned above, the weak nuclear capture rates 
strongly depend on spin factors within the total ppp 
molecular spin function and can be expressed as 

^OM = 27Qjj^(|Ag -I- jArp), 

^PM = 27pjjj(jAg -I- |Arp). 

The molecular overlap factors are 27 qj^ = 1.009 ± 0.001 
and 27 pjjj = 1.143 ±0.001 [21]. Based on these equations 
the capture rates of the molecular states can be calcu¬ 
lated using the MuCap result for Ag and the theoretical 
value for the smaller rate Arp as input (see Table |^. 

In the presence of Z > 1 chemical impurities, the muon 
can form a bound Z/r state instead of a pp atom. The 
factor / in Fig. characterizes the initial population of 
Z/i atoms, which arises from two pathways. First, at 
the time of the muon stop, Z > 1 elements are energet¬ 
ically favored over hydrogen by Coulomb capture. Sec¬ 
ond, during the pp deexcitation cascade, prompt transfer 
to higher Z elements can occur. The size of / scales lin¬ 
early with the relative atomic concentration of the 
impurity. 

Muons will also transfer from the singlet pp state to the 
energetically favorable Z/r state in collisional processes. 
Transfer from the molecular states to the Z^ state is not 
possible because the charged {ppp)^ molecule is repelled 
by the Z nucleus. The effective transfer rate to the im¬ 
purity, Ap 2 , is expressed as 

^pz = j (5) 

where is the normalized transfer rate. Excited 
Z/r states are created by such transfers, and observable 
muonic X-rays are emitted during the subsequent deexci¬ 
tation cascade. The rate A^ of subsequent muon capture 
on the nucleus increases roughly proportional to Z® (The 
more realistic Primakoff formula is discussed in mr 
Table ^ shows that the capture rates for typical impurity 
elements (nitrogen, oxygen, argon) are all much larger 
than the pp singlet capture rate Ag. 

The natural abundance of deuterium in hydrogen gen¬ 
erally causes an additional loss channel due to the forma¬ 
tion of dp atoms [271 ng and pdp molecules |g . For the 
presented measurement, a cryogenic distillation column 
was used to isotopically purify the hydrogen achieving 
a final deuterium concentration of less than lOppb |1]. 
At this level, the deuterium loss channel is completely 
negligible. 

The muon can decay from any of the states in Fig. [2 at 
a rate close to the free muon decay rate, A_|_ [HI [13]. The 
actual decay rates are slightly reduced with respect to A_|_ 
by the Huff factor h [TS] which accounts for bound-state 
corrections arising from Coulomb and relativistic effects. 
We neglect the Huff factor in the pp system in the fol¬ 
lowing equations, as it is calculated to reduce A^ by only 
26 ppm [23J[30|; in the final evaluation of Ag, Eq. 0 , we 
will explicitly include this reduction. For the argon sys¬ 
tem we use h = 0.985 ± 0.003, based on extended-model 
calculations [HIID] which include a more accurate treat¬ 
ment of finite nuclear size effects. 
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B. Kinetic Equations 

The kinetics scheme in Fig. corresponds to a sys¬ 
tem of coupled linear differential equations for the time- 
dependent populations npm(i) and n^p{t) 

of the pfi, ppp, and Z/i states. It is convenient to first 
define the total muon disappearance rate from each state: 


energy spectrum deviates from a pure Michel spectrum 
due to Coulomb effects da ISO]. 

The distribution of muon capture products (i.e., recoil 
nuclei or neutrons) versus time is 

Nc{t) = ec[Agnp^(t) -k 

^OM^om(^) + ^PM^pm(^)] + 

• (10) 


l"pAi = + ^pZ + Ag -I- A 

l^om — ^op -^OM 

Tpin = -|- Ap]^ 

Tzp = + ^z- 


ppfl 


Here Cc and account for the different efficiencies in de¬ 
tecting reaction products from capture on protons versus 
capture on nuclei with atomic number Z. 

The time distribution of X-rays from muon transfer is 

N^{t) = e^P^Ap^rip^it) , ( 11 ) 


These rates are also the eigenvalues of the system. The 
populations of the muonic states are then described by 
the following differential equations: 


n'ppit) = -Tp^ Uppit) 

n'omit) = Ki npp{t) - r, 


om ‘'om 


it) 


^pm(^) -Apf ^pp(^) T Aqp Fpjjj ripjjj(t) 


(7) 


Zu(^) — Ap2 n At) — r 


Zp ’^Zp 


At)- 


where Px is the probability for X-ray emission per trans¬ 
fer and tx is the X-ray detection efficiency. The observ¬ 
ables in Eqs. ( |10[ ) and ( |11[ ) provide the primary tools 
for experimentalists in disentangling the rich physics of 
muon-induced processes in hydrogen. 


C. Present Experimental Knowledge of the 
Molecular Formation Rate A„„,, 


The initial conditions at t = 0 are ^^^(0) = 1 — /, 
nz^(O) = / and nQjjj(O) = npjjj(O) = 0. If the rates are 
time independent, Eq. Q defines a system of differential 
equations with constant coefficients which has straight¬ 
forward but lengthy analytical solutions ni(t) (given in 
the appendix). Formally they can be written in terms of 
the eigenvalues F^ (see Eqs. (©): 

’ (^) 
cx. 

where i,a G {pp, om, pm, Z/i). Usually this is a good 
approximation, but as will be explained later there are 
cases where epithermal pp atoms are depopulated at 
energy-dependent rates in the period before they have 
fully thermalized (c.f. [sjini]). Of particular relevance 
for the present work, a muonic X-ray measurement m 
observed that the muon transfer rate ApAr(t) increased 
until it reached its constant value for the thermalized 
atom. Because Xpz{t) is time dependent, Eq. (0 must 
be numerically integrated. 

From the muonic state populations n^{t) we can derive 
the time distributions of various experimentally observed 
final-state muon-disappearance products. The distribu¬ 
tion of decay electrons is given by 


Ne{t) = A+jce [Uppit) 

e' h n- 




.(i)] 




(9) 


where Cg is the detection efficiency for electrons produced 
by muon decay from the hydrogen bound states. Depend¬ 
ing on the experimental setup, the detection efficiency e^. 
in higher-Z atoms can be different because the electron 


The basic experimental technique for measuring the 
molecular formation rate Xppp is to introduce an impu¬ 
rity to the pure hydrogen target. Though it might seem 
counterintuitive, adding this complication is helpful be¬ 
cause it opens a competing channel to molecular ppp for¬ 
mation. Since muon transfer to the impurity only pro¬ 
ceeds from the pp atom, the Z/r population follows the 
time evolution of the pp population which feeds it and the 
electron distribution described in Eq. ([^ depends mainly 
on Appp, Ap 2 , and A^. By adding the proper amount of 
a well-chosen impurity, the terms in Eq. 0 will differ in 
their time dependencies and relative sizes such that indi¬ 
vidual rates can be disentangled via a fit to the observed 
electron time spectrum. 

An early measurement of Xppp used an LH 2 target with 
deuterium admixtures |33j . In this case, muons transfer 
from pp to dp and deuterium essentially plays the role of 
the impurity Z in Fig. [^and Eq. 0- The formation of 
dp atoms can lead to muon-catalyzed fusions which emit 
gammas. Observation of the gamma yields for various 
deuterium concentrations thus enabled a determination 

of Appp. 

Other experiments employed a similar strategy. 
Ref. 1311 measured the muonic X-rays emitted follow¬ 
ing transfer to Ne. Experiment [35] simultaneously ob¬ 
served the time distribution of /iXe deexcitation X-rays 
and muon decay electrons. The first measurable deter¬ 
mined Fpp, while the second enabled independent extrac¬ 
tion of Appp and Ap^ at a single impurity concentration. 

The most recent experiment |36j used a very different 
experimental setup consisting of a layer of solid hydro¬ 
gen with various tritium admixtures. Fusion products 
were observed, and muon transfer to tritium changed the 
disappearance rate of the pp state according to the first 
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of Eqs. Q. Conceptually the experiment was therefore 
quite similar to m- 

Figure plots the relevant experimental and theoret¬ 
ical determinations of including that presented in 

this paper. The experimental data are not completely 
consistent. The higher A^p^ value measured in the solid- 
target experiment could originate from comparatively 
slower thermalization of the pfj, atoms via elastic colli¬ 
sions with the solid hydrogen lattice [32 . Review |2 
excluded the solid-hydrogen result to obtain the experi¬ 
mental world average App^ = (2.3 ±0.2) x 10® s“^, where 
the uncertainty has been inflated to account for the in¬ 
consistencies among the contributing measurements. 



FIG. 2. (Color online) Comparison of a theoretical calcu¬ 
lation 12 and experimental measurements of the molecular 
formation rate Appp. Red squares denote liquid-hydrogen tar¬ 
gets (Bleser et al. [33] ^ Conforto et al. mi the green cross 
denotes a solid-target measurement (Mulhauser et al. |36p. 
and the blue circles denote two measurements in a gaseous 
hydrogen environment (Bystritsky et al. |35| and this paper). 
The shaded region corresponds to an updated world average of 
the experimental results, excluding the outlying solid-target 
data point. 


case in which only the pp state is populated. In that case 
the electron time distribution Eq. (|^ simplifies to 


Ne{t) (X ( 12 ) 


and Ag can be determined from the difference A_ — A_(_. 

In reality, experiments must always account for ef¬ 
fects arising from the existence of muonic molecules. 
The lifetime method was pioneered by an experiment at 
Saclay m which used an LH 2 target {(j) = 1); the full ki¬ 
netics of Eq. Q therefore needed to be considered, and 
this led to significant uncertainty in the interpretation 
of the experiment’s results. The MuCap experiment jj] 
used a low-density hydrogen target (^=0.01) in order to 
more closely approach the ideal case of a purely pp, sys¬ 
tem. In the following we analyze the impact of muon 
chemistry on the lifetime method only; the reader is re¬ 
ferred to review |2 for a more comprehensive treatment 
of muon capture experiments in hydrogen. 




t ips) t (ps) 

FIG. 3. (Color online) Calculated muonic-state populations 
for (left) the hydrogen density in the MuCap experiment, 
</>=0.01, and (right) LH 2 , 0=1. In MuCap 97% of all cap¬ 
tures proceeded from the pp singlet state, while in LH 2 cap¬ 
ture takes place predominantly from ppp molecules. 


D. Impact of Molecular Effects on Muon Capture 
Experiments 


Muon capture experiments determine Ag either by 
measuring the rate of neutron emission according to 
Eq. (101 (“neutron method”) or by inferring the muon 


disappearance rate in hydrogen, A_, from the time dis¬ 
tribution of electrons, Eq. ^ (“lifetime method”). While 
the neutron method does not require high statistics, its 
precision is fundamentally limited by the fact that the 
neutron detection efficiency Cc must be known to a level 
that is difficult to achieve in practice. Conversely, the 
lifetime method requires high statistics but absolute de¬ 
tection efficiency is not a factor. The basic idea of the life¬ 
time method can be illustrated by considering the ideal 


Figure shows the time distributions of pp and ppp 
populations in the hydrogen targets used in the Mu¬ 
Cap |3] and Saclay m experiments. At the lower target 
density used in MuCap, muons remain predominantly in 
the singlet pp state over the course of the typical mea¬ 
surement period of 15 microseconds. There is neverthe¬ 
less non-negligible formation of (ppp)^^ molecules, and 
therefore good knowledge of the rate App^ of the process is 
necessary for correct interpretation of the experiment. In 
contrast, in the LH 2 target used in the Saclay experiment 
the muon quickly populates the {ppp)^^ state, within 
1 ps, and the subsequent depopulation of the (ppp)^^ 
state to the {ppp)^^ state at rate A^p is the crucial ele¬ 
ment to interpreting the experiment. 
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III. EXPERIMENTAL METHOD 


B. Measurement of Aj,p^ in Argon-doped Hydrogen 


A. MuCap Apparatus 


The MuCap detector (Fig. will be described here 
only in brief; greater detail is available in [iiiiiMiEn]. 
The experiment was located at the TrifS secondary muon 
beamline of the 590 MeV proton cyclotron at the Paul 
Scherrer Institute. Low-energy muons (34 MeV/c) passed 
through a scintillator counter (/iSC) and a wire-chamber 
plane (/iPC) before coming to a stop inside a 10-bar hy¬ 
drogen time projection chamber (TPC). 

The ^SC provided the start signal for the muon life¬ 
time measurement, and the fiSC and ^PC together pro¬ 
vided efficient pileup rejection which enabled selection of 
events in which only a single muon was present in the 
TPC. The TPC |301 01] provided tracking of incoming 
muons and clear identification of each muon’s stopping 
location by detecting the large peak in energy deposition 
at the end of the muon’s Bragg curve. The trajecto¬ 
ries of outgoing decay electrons were reconstructed by 
two concentric multi-wire proportional chambers (ePCl 
& ePC2), while a scintillator barrel (eSC) provided the 
stop time for the lifetime measurement. 



FIG. 4. Simplified cross-sectional view of the MuCap detector 
setup. Neutron detectors not shown. The main components 
are described in the text.(Figure reproduced from |10|.l 


Fiducial cuts can be applied to the TPC data to se¬ 
lect muons that stopped in the hydrogen gas, far away 
from any rate-distorting Z > 1 materials. The three- 
dimensional electron tracking makes it possible to corre¬ 
late a decay electron with the stopping point of its parent 
muon, thereby increasing the signal-to-background ratio 
(see [^133] for details). 


For this measurement we introduced argon to the oth¬ 
erwise ultra pure hydrogen gas, which was at density 
(j) = 0.0115 ± 0.0001. The atomic concentration of argon 
was = 19.6 ± 1.1 ppm, as measured both volumet- 
rically during the initial filling and by gas chromatogra¬ 
phy at the end of the measurement. The two concen¬ 
tration measurements were consistent, but we have con¬ 
servatively expanded the uncertainty to cover the uncer¬ 
tainties of both. The gas density has been derived from 
the temperature and pressure, which were continuously 
monitored. 


In principle, two time distributions, Nc{t) (Eq. @), 
and Ne{t) (Eq. @), are experimentally observable. The 
former can be measured either by using the TPC to de¬ 
tect nuclear recoil signals from /r -I- Ar —> Cl* -I- v capture 
events, or by using liquid scintillators to detect neutrons 
emitted by the excited final-state nucleus. There are two 
disadvantages to measuring Nc{t). First , the spectrum 
determines F^^ (c.f. in Eqs. (All) and therefore 

only the sum of the two transfer rates A^^^ and not 

the individual rates themselves, and A^^,. is not known 
with sufficient precision to enable A^p^ to be extracted 
independently. Second, there are significant systematic 
uncertainties relating to spatial pileup of TPC signals 
from the stopping muon and the capture recoil, and to 
uncertainties in the neutron time of flight. 


The MuCap experiment was designed to detect decay 
electrons, so we used a high-statistics sample of Ne{t) to 
extract Xpp ^. If a muon decays it cannot undergo nuclear 
capture, eliminating the possibility of distortions in muon 
stop identification due to additional energy deposit from 
capture recoils. Consequently, the analysis and system¬ 
atic uncertainties were very similar to those developed 
for the earlier lifetime experiment measuring Ag 0 ]. 

The decay-electron analysis works as follows. With the 
judicious choice of argon concentration c^j. = 0(20 ppm), 
the disappearance rates Fp^ and F^^.^ in Eq. ([^ are suf¬ 
ficiently different as to allow them to be unambiguously 
extracted from a fit to the corresponding decay-electron 
time spectrum. The argon capture rate A^,. [14] is three 
times higher than the muon decay rate and therefore 
transferred muons disappear quickly. Under our condi¬ 
tions, the contributions of App^ and Ap^,. to the total 
(p/i)g disappearance rate Fp^ were 4% and 8%, respec¬ 
tively. As above, the eigenvalue Fp^ alone would only de¬ 
termine the sum of two unknowns, App^ and Ap^,.. How¬ 
ever, both rates enter into the coefficients c“ in Eq. ([^ in 
independent combinations, as can be seen from the full 
solutions in the appendix. A combined fit can therefore 
simultaneously determine App^, Ap^,., and, as a byprod¬ 
uct, A^,., without any need for absolute normalization. 
To address concerns about the uniqueness and stability of 
this multi-parameter fit to a single distribution, we per¬ 
formed extensive pseudo-data Monte Carlo studies of the 
full kinetics equations; good convergence was observed. 
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IV. ANALYSIS AND RESULTS . ^ ^ , n . 

TABLE II. Experiment specific parameters used in the fat of 

Eq. ( |13[ l to the data. See text for detaiis on their evaiuation. 

A. Data Analysis 


A total of 7.2 X 10® fully reconstructed muon decay 
events were used in the present analysis. These events 
were selected via application of our standard cuts, de¬ 
scribed in Ref. [1]. Each event was required to involve 
a pileup-free muon stop in the TPC fiducial volume, 
AxxAyxAz = 10.4x8.0x20.4 cm®. The decay-electron 
trajectories were reconstructed from spatial and tempo¬ 
ral coincidences among the two cylindrical wire chambers 
and the two layers of plastic scintillators. Once the set 
of good events had been selected, the time differences 
between the fast signals of the electron scintillator eSC 
and the muon beam scintillator fj,SC were histogrammed 
and the resulting decay time spectrum was fitted with 
the function 


N{t) = A [np^(t) + np^(t) + n^^(t) + e h + B 

(13) 

using the MINOS package. This fit function is identical 
to Eq. ([^ apart from the introduction of a flat back¬ 
ground term B. The relative efficiency e is defined as 
e = e'Jee- 

To accommodate the time dependence of ApAr 
in a nearly model-independent way, this rate was 
parametrized in the form 


= ApAr(l - 0.6 


(14) 


where a and /3 were extracted from Fig.l in |17j . The pa¬ 
rameter /3 characterizes py thermalization and was scaled 
down by 1.5 from the value in as that experiment 
used a 15-bar target whereas MuCap used a 10-bar tar¬ 
get. The scaling of a with pressure depends on the ini¬ 
tial population of hot pfi atoms after the muonic cascade, 
which, according to theory m. should increase by ~10% 
with a pressure increase from 10 to 15 bar. We did not 
change the value of a extracted from m, but we as¬ 
signed it a conservative 50% uncertainty. The final fit 
method used numerical integration with the values listed 
in Table The analytical solution (A21 was used for 
cross checks. 


The fitting procedure using Eq. (13) requires a tim¬ 


ing calibration to assert that the muon arrival time is at 
t = 0. For that, the rising edge of the histogrammed 
differences of the pSC and the sixteen eSC subdetectors 
were fitted individually. This determined timing calibra¬ 
tion offsets for each eSC detector with a precision of 2 ns. 
The sixteen offsets were then applied to their correspond¬ 
ing spectrum before the sum of all time distributions was 
fit with Eq. (131. 


The fit was performed over the range [0.12/xs, 20^s]. 
Five quantities were treated as free parameters: A^^^, 
A^,., the normalization A, and the background 
term B. All other parameters were fixed in the fit to 
the values in Table |T] and, for experiment specific pa¬ 
rameters, according to the values given in Table |TT] The 


Parameter Value 


CAr 19.6 ±1.1 ppm 

(j) 0.0115 ±0.0001 

/ 5 ± 1 X 10~^ 

e 0.996 ± 0.003 

Cq 57 ± 57 ppb 

Cn 115 ± 115 ppb 

a 0.25 ±0.12 

P 1.0±0.2 X 10'^ 


initial p,Ar formation fraction / = (5 ± 1) x 10“'^ is the 
sum of two components, fc and fe- The atomic capture 
ratio for argon relative to hydrogen has been measured 
to be /c = (9.5±1.0)cAr = (1.87±0.20) x An 

additional initial population /e = (1.66 ± 0.34)/c from 
excited-p/i-state transfer has been observed in a target 
at 15-bar pressure |13]. We account for this by using 
/e = (3.1 ± 0.9) X 10“"*^, in which the uncertainty has 
been conservatively enlarged to accommodate the possi¬ 
bility of a pressure dependence. 

The energy spectra of decay electrons emitted from pp 
and Arp atoms are different, which leads to a difference 
in the corresponding detection efficiencies. We used the 
energy spectrum calculated in gj] and folded it together 
with the energy-dependent detector efficiency obtained 
from a full Geant4 simulation. The resulting relative ef¬ 
ficiency, e = 0.996 ± 0.003, shows that the thin layers of 
the MuCap electron detectors are not very sensitive to 
spectral differences at higher energies. 

After the fit, small corrections were applied to the fit¬ 
ted rates to account for the presence of chemical impuri¬ 
ties oxygen and nitrogen, with atomic concentrations Cq 
and CjY, respectively. This procedure is discussed in the 
next section. 


B. Results and Systematic Uncertainties 


The fit to the data is plotted in Fig. The upper 
panel shows the decay electron time spectrum alongside 
the time distributions of the parent muon populations 
iT'Pfj.it), ^-Ar^)^)’ '^om(^)> and Tip^it) determined by the 
fit. The lower panel displays the residuals, i.e., the differ¬ 
ences between the data and the fit function normalized by 
the uncertainty of each data point. The good agreement 
between the data and the fit function is demonstrated by 
the reduced y^/DOF = 0.983 ± 0.064. 


Table III presents the fit results for the three rates 
App^, ^pAn and A^^- The table also lists systematic 
corrections A and the systematic uncertainties S resulting 
from a Ala variation of the fixed parameters listed in 
Tables U and HH 

The fit did not explicitly model effects from the accu- 
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FIG. 5. (Color online.) Upper Panel: Fit to the decay elec¬ 
tron time spectrum using Eq. ( |13[ |. The data is shown as the 
black line. The colored curves depict the time-dependent con¬ 
tributions from the kinetic states; the black dashed line is the 
fitted sum. Lower Panel: The normalized residuals between 
the data and the fit function, {Ni — N{ti))/ai, indicate good 
agreement in the fitted range 0.12-20 /rs. 


TABLE III. Fit results for A^^r) S'lid A^r, as well as 

their associated systematic corrections (A) and uncertain¬ 
ties (5). The final error on each rate is the quadrature sum 
of the contributing uncertainties. 




[s-^1 

ApAr [S“^] 

Aa, [lO^s-'l 

Fit 

22 996 

±647 

43 799 ±151 

13 023 ±147 

Systematic 

A 

5 

A S 

A S 

Timing calibration 


39 

13 

35 

Efficiency e 


37 

13 

34 

Huff factor h 


45 

15 

27 

f 


46 

12 

39 

As, Aqmi ApM 


13 

16 

5 

-^op 


31 

3 

2 

Apf 


9 



Epithermal 


329 

85 

278 

H 2 O and N 2 

116 

116 

-15 15 


Final result 

23112 

±741 

43 784 ±177 

13 023 ±322 


mulation of nitrogen and oxygen in the hydrogen due to 
outgassing from the TPC vessel. Instead, a correction A 
was applied to the fitted values of both A^p^ and Ap^^.. 
During its main run MuCap achieved hydrogen chemical 
purity levels of better than 10 ppb, but during the argon- 
doped measurement the TPC was disconnected from the 
hydrogen circulation and purification system |45j . After 


six days, atomic concentrations of Cq=115 ppb of oxy¬ 
gen (in the form of water vapor) and Cf^=230 ppb of 
nitrogen were observed using a humidity sensor and gas 
chromatography, respectively. Due to the higher muon 
transfer and capture rates for oxygen compared to ni¬ 
trogen, transfer to oxygen is the dominant effect need¬ 
ing to be taken into account in the correction to the 
measured rates. A series of pseudo data were gener¬ 
ated based on the kinetics in Eq. Q, with transfer to 
and capture on chemical impurities included. MuCap 
had previously measured these rates independently us¬ 
ing impurity-doped hydrogen mixtures. Our result for 
ApN agreed with previous measurements, but our results 
for Apo (measured via water doping) were nearly two 
times higher than the value quoted in Table |Tj For inter¬ 
nal consistency we used the transfer rates measured by 
MuCap in our simulation. The pseudo data were then 


fitted with Eq. (131 to extract the shifts in the rates 
Appp, Ap^j., and as a function of the oxygen con¬ 
centration Cq. As the exact time dependence of the im¬ 
purity buildup was unknown, conservative estimates of 
Cq = 57 ± 57 ppb and Cj,j = 115 ± 115 ppb were used to 
cover all possible accumulation scenarios. The impurity- 
related corrections to Ap^^ and Ap^^j. were determined to 
be App^ = 116 ± 116 s“^ and Ap^^. = —15 ± 15 s“^. 

The final results for the fitted rates after applying the 
impurity-related corrections and summing all systematic 
uncertainties (Table III I are 


= 2.311 ± 0.074 X 10'‘ s'^ 

ApAr = 4.378 ± 0.018 x lO^^ s'^ 

Aa^ = 1.302 ±0.032 X 10® s'^ . 
From these one can deduce the normalized rates 


(15) 


= 2.01 ± 0.07 X 10® s-i 

ApAr = 1-94 ±0.11 X 10“ s-i (16) 

using Eqs. (§ and ([^, respectively. 

The normalized correlations among the five free fit pa¬ 


rameters are presented in Table IV These correlations 
are incorporated into the uncertainties on the final re¬ 
sults. 


TABLE IV. Normalized correlation coefficients of the free pa¬ 
rameters in the fit to the decay electron time spectrum. 


Rates 

■^ppp -^pAr 

Aap a 

■^pAr 

0.9548 


Aaf 

-0.8021 -0.9011 


A 

0.0495 0.0269 

0.0234 

B 

-0.6603 -0.5479 0.4189 -0.1082 


Our result for App^ is about Icr larger than the value 
we obtained in [3] due to the more refined analysis in 
this paper and the correction of a numerical error in the 
fitting code. As regards the transfer rate to argon ApA^, 
there is a wide spread of experimental results obtained 
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with different methods and target conditions, clustered 
around 1.4 x 10^^ s“^, 3.6 x 10^^ s“^, and 9 x 10^^ s“^, 
as discussed in m- Our value = 1.94 ± 0.11 x 
10^^ s“^ is close to the most recently published value, 
1.63 ± 0.09 X 10^^ s“^ [13, albeit 2.2cr higher. Note that 
the uncertainty in the argon concentration only enters 
into the extraction of the normalized rate while in 

the fit to determine A^^^ effective rates are being used 
which are independent of CAr- Our result for the muon’s 
nuclear capture rate on argon, A^^., agrees well with the 
values in the literature, 1.20 ± 0.08 x 10® s“^ |3B] and 
1.41 ±0.11 X 10® s"^ @7] . 


C. Consistency Checks 


The fit start time was varied to check for any distor¬ 
tions or physical effects not accounted for by the fit func¬ 
tion. Figurej^shows the progressions of the fitted rates as 
the fit start time was increased in steps from its standard 
value of 0.12 ^s. The red lines denote the ±1 ct variation 
allowed because of the set-subset statistics involved in 
this procedure. Each rate is statistically self-consistent 
across the fit start time scan. 

In the fit to the data using Eqs. and ( [I^ , the cap¬ 
ture rate Ag is required as an input to extract A^p^, while 
the latter is itself used in the determination of Ag. This 
interdependency is not a problem because all fitted rates 
in Eq. (151 depend only very weakly on the hydrogen 


capture rates, as quantified in Table |III| We explicitly 
iterated the procedure (obtaining fit results with Ag as 
input, using the results to correct Ag, repeating with the 
adjusted Ag) to arrive at a stable, self-consistent solu¬ 
tion, and we found that the results for App^, Ap^,,, and 
■^Ar were changed by less than one tenth of their uncer¬ 
tainties. 

Lastly, the reproducibility of the fit was tested by gen¬ 
erating 10"^ pseudo-data histograms using the final fit pa¬ 
rameters in Eqs. (15), and fitting each pseudo experiment 


in the same manner as the real data. The fits consistently 
yielded the input values, and the simulated data repro¬ 
duced the same fit uncertainties listed in Table [ 


V. RELEVANCE TO THE INTERPRETATION 
OF THE MUCAP EXPERIMENT 


In Section the influence of the molecular rates App^ 
and Aqp on muon kinetics in hydrogen was described. The 
MuCap experiment measured the effective muon disap¬ 
pearance rate A_ in low-density ultra-pure hydrogen by 
fitting the observed decay electron time distribution with 
a three-parameter function, f{t) = ± B. Taking 

ppfi formation into account, the disappearance rate can 
be expressed as 


(17) 


Start time scan App^ 



Start time scan 



44200 
i 44100 
44000 
43900 
43800 
43700 
43600 
43500 
43400 



FIG. 6. (Color online) Results from a scan over the fit start 
time. Each plot shows the fit results (blue points) for a par¬ 
ticular rate as the start time is varied. The variation of each 
rate is consistent with the expectations from the one-sigma 
statistically allowed set-subset deviation (red line). 


Here AA„,, = —12.3 s ^ is a calculable bound-state modi- 
fication to the muon decay rate in the pp system |29II30] . 
and AAppp is a modification to Ag accounting for the 
small population of muonic molecules and the fact they 
have unique capture rates. In the following we summarize 
the derivation of AApp^, based on our improved measure¬ 
ment of Appp at conditions nearly identical to those of the 
main MuCap experiment. 


The derivation is based on high-statistics simulations 
of the full kinetics described by Eqs. (A1). Since the Mu¬ 
Cap measurement of Ag was performed using pure hydro¬ 
gen gas, for the simulations the Z channel was used to 
model the small amounts (few ppb) of oxygen and nitro¬ 
gen impurities that were observed to have outgassed from 
the hydrogen vessel’s walls. The relevant input parame¬ 
ters for the simulation were those in Tables in and ITU An 
accidental background was added to make the signal-to- 
background level commensurate with that in the MuCap 
data. Time distributions of 10^^ decay electrons were 


A- — Ag. ± AAp^ + As ± AApp^ . 
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generated for two different cases: = 0 and 

The previous MuCap analysis |1] was performed using the 
preliminary value A^'^^'^p = 1.94 ± 0.06 x 10® s“^; here 
we update the analysis using our new result in Eq (211. 
To determine the effect on the MuCap result for Ag, we 
fit the simulated time distributions with the same three- 
parameter function used to fit the data. The relevant 
correction is then obtained via 


= A_ (A“f ^P) - A_ (Appp = 0) , (18) 

where the A_ values are obtained from fits to the two 
simulated data sets generated using different A^^^ val¬ 
ues. The uncertainty in AA^^^ is estimated in a similar 
manner, by generating pseudo data while varying the pa¬ 
rameters entering the kinetic equations by ±lcr individu¬ 
ally. The resulting fit determines the final correction for 
the MuCap experiment to be AAp^^ = —18.4 ± 1.9 s“^, 
which is smaller than the correction in a by 0.7 s-i. 
Thus the updated value of A“p®^**'P induces a small shift 
of the singlet pfj, capture rate measured by MuCap from 
Ag = 714.9 ± 5.4stat ± S.lgyst obtained in a to 

Ag = 715.6 ± 5.4stat ± 5.1,yst . (19) 

The value of the pseudoscalar coupling constant, gp = 
8.06 ± 0.55 extracted in a) is correspondingly changed 
by —0.045, i.e. by only 8% of its uncertainty. 

From our simulations we can determine the depen¬ 
dence of AApp^ on molecular parameters, 

= -18.4 [1 -I- a(App^ - A“p®^*"P) -I- - A°p)] , 

( 20 ) 

where A^p is given in Table hi a = 4.7 x 10 and b = 
2.9 X 10“®. Using the new measurement presented in 
this paper, the total uncertainty in the MuCap capture 
rate Ag due to ppg formation is less than 2 s“^ and is 
dominated by A^p, while App^ contributes only 0.6 s“^. 


Our result for the ppp, formation rate, 

A“"pC^P = 2.01 ± 0.06stat ± 0.03,yst x 10® s-\ (21) 

is 2.5 times more precise than previous measurements, 
which were performed under a variety of different ex¬ 
perimental conditions and whose results disagreed be¬ 
yond their uncertainties. To obtain a new world average 
we used the procedures for averaging and inflating un¬ 
certainties advocated by the Particle Data Group m- 
we have included only the gas- and liquid-target experi¬ 
ments [55H55] . choosing to omit the lone solid-target ex¬ 
periment m because of possible solid-state effects which 
are not well understood. The updated experiment world 
average then becomes 

A^^« =2.10 ±0.11 X 10® s-i . (22) 

The rate App^ was a necessary input to the MuCap 
experiment’s recent precision determination of the nu¬ 
clear capture rate on the proton, Ag |1]. MuCap was 
designed so that the majority of muons underwent cap¬ 
ture in muonic pp, atoms, and formation of ppp molecules 
changed the observed capture rate by only 2.5%. How¬ 
ever, given the inconsistency between existing App^ re¬ 
sults it was difficult to confidently estimate the uncer¬ 
tainty on the correction to Ag for ppp effects. Our new re¬ 
sult for App^, obtained at the same hydrogen density and 
temperature as in the main MuCap experiment, leads to 
a well-defined correction to Ag, and the corresponding 
contribution to the total error is now minor. The value 
of presented here differs only slightly from the 

value used in [3], and consequently the updated values 
for Ag and the pseudoscalar coupling gp agree to better 
than O.Ict with the values in that publication. 


VI. SUMMARY 

The time spectrum of electrons emitted by the decay of 
muons stopped in argon-doped hydrogen were measured 
with the MuCap detector, for the purpose of determin¬ 
ing the formation rate App^ of ppp muonic molecules. 
The TPC enabled selection of muons that stopped in 
the hydrogen, away from high-Z materials, and the elec¬ 
tron tracker provided Stt solid-angle coverage and enabled 
vertex matching with muon stops. We developed a de¬ 
tailed physics model to describe the time evolution of the 
atomic and molecular muonic states contributing to the 
decay electron spectrum, taking into account the energy 
dependence of the muon transfer rate Ap^,. from hydro¬ 
gen to argon. We extracted App^, Ap^,., and the muon 
capture rate in argon, A^,., from a single fit to the decay 
electron time spectrum. Our results for Ap^,. and A^,. 
agree with those from previous dedicated experiments. 
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Appendix A: Solutions to the muon kinetics 
equations 

The differential equations in Eq. 0 can be solved by 
determining the eigenvalues and eigenvectors of the sys¬ 
tem. The time-dependent populations of the four muonic 
states are given by 


it) = (!-/)• 


PM 




A, 


of 


r -r 

om pfi 




At) = 


1 -/ 


pmv''/ -p -p \ T'' I”' 

^ om ^ pp \ ^ om ^ pm 


Aof^op _ g-rp„t^ r^^Apf ApfA^p ^ 


r - r 

pp pm 


(Al) 




W = (1-/) 


^pZ 

r — r 

^ Zfj. ^ pp 


(e —e ^2p*)-|-/-e . 


One simplistic but heuristic approximation is to neglect 
the small parameters Ag and Apj^ in the disappearance 
rates in Eqs. ([^ and the initial Z/r population /, and 
assume that T^ is large compared to all other eigen¬ 
values in Eq. (1^. In this limit the observable electron 
distribution, Eq. attains the simple form 


Ne{t) oc e-^+* 


1 + 

A™,, 


(A2) 


which elucidates our strategy of determining A^^^ in a 
single fit. 
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